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Abstract 



In this article we present a continuous time model for natural gas and crude oil 
future prices. Its main feature is the possibility to link both energies in the long 
term and in the short term. For each energy, the future returns are represented as 
the sum of volatility functions driven by motions. Under the risk neutral probabil- 
ity, the motions of both energies are correlated Brownian motions while under the 
historical probability, they are cointegrated by a Vectorial Error Correction Model. 



Our approach is equivalent to defining the market price of risk. This model is free of 



arbitrage: thus, it can be used for risk management as well for option pricing issues. 
Calibration on European market data and numerical simulations illustrate well its 
behavior. 

Key words: future prices, natural gas, crude oil, cointegration, Vectorial Error 
Correction Model, arbitrage free modelling. 
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Introduction 



The need to model future prices of gas and oil simultaneously to optimize 
energy portfolios is now very present. Indeed an energy portfolio could be 
invested into several energy markets which could interfere with each other. 
For example, a firm which detains a supply contract will need to model gas 
and oil evolution for an optimal risk management. 



With the energy market liberalisation over the last decade, one may think that 
gas and oil prices are decoupled, but actually several statistical studies tend 
to prove that prices are cointegrated. For statistical evidences of cointegration 
and economic explanations, we refer for instance to the articles by Panagiotidis 
and Rutledge 15] and Asche et al. [l| for the UK and European markets, and 
to the article by Bachmeier and Griffin ^ for the US market. The dependence 
between gas and oil prices could be economically explained with gas long term 
contracts, which still represent the majority of supply in European gas and 
whose prices are indexed on oil and oil products prices. This indexation creates 
a structural link between prices of both energies. 



In the previously cited references, one proposes econometric models (in dis- 
crete time) for gas/oil prices, which are coherent with market data and which 
accounts well for the interdependence of prices. They are useful for some risk 
management purposes, such as Value at Risk measurements. However, as soon 
as we have to consider energy contracts and related pricing/hedging issues, dif- 
ferent models emerge in order to be consistent with the arbitrage free theory: 
they are such that forward contracts are martingales under risk neutral prob- 
abilities (see Musiela and Rutkowski [l3| ) . In the following and to simplify our 
presentation, we identify future prices (given by quotations data) and forward 
prices (given by models), which is correct if interest rates are deterministic for 
instance. Usual factor models for the spot and forward contracts on a given 
energy are written as 

dF{t,T) 



F{t,T) 

where 



a(t,T)mB t 



• F(t, T) is the forward contract quoted in t and delivered in T, 

• cr(t, T) = (o"i(t, T), o"2 (t, T), . . . , cr„(t, T)) is a row vector of normalized volatil- 
ity functions of forward returns, 

• n is the number of risk factors identified through a PCA (Principal Com- 
ponent Analysis) of forward returns, 

• B t = (Bj, B^, . . . , B™)* are independant Brownians motions under risk neu- 
tral probabilities (here * stands for the transposition), 

• S is a n x n matrix, equal to the square root of a variance-covariance matrix. 
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See the works by Geman Clewlow and Strickland \^ among others. 

The volatility functions describe the shifting, the twisting and the bending 
(for Clewlow and Strickland 0). In Brooks [3], they describe the level, the 
slope, the curvature (LSC model). R. Brooks uses this model for natural gas 
contracts and gives an explicit form for these functions : 

a\{T - t) = 1, 

a e 2 (T-t) =e~¥, 
o-UT -t) = 1 —t e ~^i for % > 3. 

The model is linear in the parameters (E and (rf )j), thus ordinary least squares 
regression is applied to estimate them. 

For crude oil contracts, analogous models could be set up, with different Brow- 
nian motions. These ones can be correlated to those of gas models. Within this 
approach, it appears that volatility functions adjust well for each energy. Nev- 
ertheless, long term dependences are poorly modeled. The variance-covariance 
matrix (related to S) induces relevant marginal distributions of the returns 
but unfortunately unrealistic joint price distribution. For example, when sim- 
ulating the model one often obtains a growing prices' scenario for gas and a 
decreasing one for crude oil (see Figure 4 in Section 3). It is fundamental to 
note that the previously decribed models (1) are written under the risk-neutral 
probability (denoted by Q in the sequel), which is suitable for pricing/hedging 
issues, while the long-term dependence (given by an econometric cointegra- 
tion analysis) holds under the historical (or physical) probability (denoted by 
P). To accommodate both features (long term dependences on the one hand; 
stochastic returns described by volatility functions on the other hand), a natu- 
ral idea consists in suitably modeling the market price of risk (Xt)t, making the 
connection between historical and risk-neutral worlds. This is the main contri- 
bution of our work. Details are given in Section 2. Thus, our model cointegrates 
gas and oil prices while being coherent with the pricing by arbitrage. In the 
following, we mainly focus on natural gas and crude oil. 

We now mention a similar approach to ours. In his PhD thesis, Steve Ohana 
aims at modeling spot and forward contracts for two cointegrated energies, 
but in discrete time (see Ohana [13]). He handles the case of US natural gas 
and crude oil market. For each energy e (e = g for g clS, 6 — C for crude), the 
model is given by (under P) 

AJ 7*' r) = e-^-O AX* + AYf 
F e (t, T) 1 1 

where is the characteristic time of the short term shock. The vector (AXf , AX£, 
AY t 9 , AY t c ) has a drift component, equal to a constant plus two terms related 
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to the past and present of the process, and a noise component expressed as 
independent GARCH processes. The dependence on present is defined by the 
non-linear link between (AVf , AXf) and the processes (Xf ) and (V t c ) and 
between (AY t 9 ,AY t c ) and the processes (Y t 9 ) and {Y t c ). 

Finally we mention very recent works, whose main motivations are to anal- 
yse the market risk premium n e (t,T) through the specification of the market 
price of risk (\t)t- We recall that the market risk premium of the energy e 
is defined by the difference between forward prices and expected future spot 
prices (under the physical probability PlPl: 



7r e (t,T) = F e (t,T) -E v (S T \J r i 



th 



where Tt stands for the information up to time t. Cartea and Williams [|| are 
dealing with the two factors long term/ short term model in the US gas market. 



Kolos and Ronn provide a statistical methodology to estimate constant 
commodity market prices of risk. Benth et al. j3j make explicit the connection 
between A and 7r e through the market players' risk preferences. All these works 
focus on an energy at a time. Our approach is rather different because we 
specify the risk market price in order to design long term dependences between 
gas and oil. Regarding the market risk premium, with a suitable choice of (Xt)t 
we propose to adjust it to at time because it is often a requirement of 
practitioners. However, it could be adjusted to any value, to generate positive 
or negative market risk premia. 

The model presented in the following (Section 2) is mainly inspired by the 



works of Brooks [4| and Ohana [14|. We use a cointegration approach and 
we design a Vectorial Error Correction Model (VECM). This type of model 
is shortly discussed in Section 1. Compared to Ohana's model, a continuous 
time setting is used, which is compatible with the arbitrage free theory. Note 
that the cointegration and VECM are now standard tools in energy markets: 
see for instance De Vany and Walls [7| for spot electricity prices, Serletis and 



Herbert 17| for US gas/oil prices. Finally, in a third part, we make numerical 



simulations to test the model. 



1 Cointegration and Vectorial Error Correction Model (VECM) 

Expressed in 1987 by Engle and Granger the cointegration concept char- 
acterizes the fact that in long term, a specific combination of non-stationary 
processes could be stationary. A time series vector {y t = [y] , . . . , yf) : t E N} 
is said cointregrated if each series (y\) is integrated with an order of integra- 
tion equal to 1, and if some linear combination of the series a -y t is stationary. 

1 note that the sign of 7r e is not related to backwardation and contango situation. 



4 



This linear combination is the long term equilibrium between each component 
of the vector y t . 



We can give a simple example based on Hamilton [10 



Ay] = AW t \ 

yl = 2y\ + AW£ V 

vl = vl = o. 



(2) 



Possible simulation of this model is given on Figure 1. If (AW/,AW t 2 ) are 
uncorrelated white noise process, {i)],yf) are cointegrated. 




Figure 1. Example of cointegrated time series. 

The graph on the right side allows us to see the scatter plot of both processes and 
the correlation between their trajectory. The size and the color of the font change 
according to the value of the correlation: the higher the correlation is, the larger the 
font is. 



The Equation (2) can be rewritten : 



^ Ayl = AW t \ 
Ay? = 2yl-y? + 2AW t 1 + AW t 2 . 



This is the form of a simple VECM (Vectorial Error Correction Model). More 
generally, if a time series vector {y t = (yj, . . . , y") : t G N} is cointegrated, the 
following relation could model the cointegration relation between the compo- 
nents of y t : 

Ay t = Uy t + EAW t . 

The matrix II represents the long term relation while the matrix E stands for 
the short term behavior. 
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We extend this notion to a continuous time framework, by writing 



dy t = Tly t dt + mW t (3) 

where (W t )t is a vector of independent Brownian motions. 

We point out that (ytjts dynamics identifies with that of a multidimensional 
Ornstein-Uhlenbeck process, with possibly a degenerate matrix EL In fact, 
several cases are possible: 

• If the rank of II is full (rk(U) = n), motions yt may be stationary processes 
and n cointegrating relations exist. 

• If the rank of II is zero (rk(U) = 0), motions y t are not cointegrated and 
not stationnary. 

• If the rank of II is r with < r < n, two (n x r) matrices a and (3 exist 
such as II = a(3* . The r linear independent columns of (3 are cointegrating 
vectors and the elements of a determine the speed of the long term return. 
Motions yt are cointegrated and non-stationary. 



For more details, see Pfaff 16 



2 The model 



To simplify our presentation, we assume that the prices of energy contracts are 
not correlated with the currencieg[f] and that interest rates are deterministic!!]- 
In the following, we choose to model gas and oil forward prices as the sum of 
volatility functions, which are driven by cointegrated motions. These motions 
depend on a deterministic function which can be adjusted to centre prices on 
the initial forward curve. The cointegrated term is related to the market price 
of risk. 



2. 1 Derivation of the model 



We first write the dynamics of the forward prices' return for both energies 
under the risk neutral probabilitjQ Q (this is the same dynamics as Brooks' 



2 To avoid extra contributions related to currency volatilities in the model. 

3 This makes future and forward prices equal. 

4 Under our assumptions on interest rates, all the forward risk-neutral probabilities 
coincide (and are equal to Q). Thus for any T, (F e (t,T)) t is a martingale under Q. 
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model): 



= *%T - t)dX t , (4) 
dX t = EdS t , 



where we set 



• F e (t, T) for the forward price quoted in t, delivering at the date T one unit 
of the energy e (g for natural gas or c for crude oil). 

• N e for the factors' number for the energy e (in practice N e = 3). This is 
the number of motions kept in the PCA during the calibration. 

• B t = (B], . . . , b^ 3+nc )* for independant Q-Brownian motions. 

• £ for a non-degenerate (N 9 + N c ) x (N 9 + N c ) matrix, used to correlate 
the returns. 

• (of (T — t))i for the set of normalized volatility functions for the energy e: 

• <ri{T-t) = l, t 

■ al(T-t) = e^, T 

• a 3 e (T-t) = ^e~¥. 

• o 9 (T — t) and cr c (T — t) for the volatilities of each energy as (N 9 + iV c )- 
dimensional row vectors defined by 

^(T-t) = K,...,< 9 ,0,...,0)(T-t), 
a\T - t) = (0, . . . ,0,<7?, . . . ,a c m ){T - t). 

(X t ) t are the so-called N 9 + A^ c motions of gas and oil: under Q, we have 
defined them as non-normalized Brownian motions. In the above definition, 
the initial value of X is arbitrary because only increments of X have an impact 
on the dynamics of forward prices (under Q) . 



In the working paper of Brooks [4], we can find a complete explanation of the 
different factors. <j\ is the volatility coefficient for the level of change in the 
forward curve. Since dF(t,T) behaves as F(t, T)o"idX i 1 when T — t tends to 
infinity, the first risk function drives the long term volatility. o"2 is the volatility 
coefficient for the slope of changes in the forward curve, and the other term 
(J3 is the volatility coefficient associated to the curvature. 

Let us go back to a model under the historical probability P. For this, we take 
advantage of the fact that (B t ) t corrected by the market price of risk (Xt)t 
becomes a P-Brownian motion, which we denote by {W t )t (see Musiela and 



Rutkowski [13j). We have 



B t = W t + f'Xsds. 
Jo 
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(Xt)t is a (N 9 + N c ) dimensional process, which we define as 



A f = Z- 1 [IlX t + r ]t } 



(5) 



where 



• II is a (N 9 + N c ) x (N 9 + N C ) matrix representing the long term link between 
the motions (X t ) t , 

• {f]t)t is a (N 9 + iV c ) -dimensional process taken as deterministic in the fol- 
lowing. This will play the role of a centring factor. 

This definition of (Xt)t is crucial in our model. In this way, A represents the 
cointegration effect between energies, by putting a long term relation between 
the motions (X t ) t : 



The term HX t cointegrates motions of gas and oil. The above equation is a 
continuous time VECM as the ones described by (3). 

2.2 Final form of the model 

To sum up, under the historical probability P, the dynamics of the forward 
prices return is: 



where W is a standard P-Brownian motion. Under (spot and forward) risk 
neutral probabilities, one has dX t = YidB t with the standard Q-Brownian 
motion B. 



2.3 Calibration of the model 

In practice, the available data are related to gas and oil future contracts, and 
not to derivatives on gas/oil (which are OTC contracts)!"^! Hence, we perform 
a statistical calibration. 

The determination of the parameters is done in five steps: 



5 Market data related to energy derivatives would be useful to better calibrate the 
volatility functions. 



dX t = TlX t dt + EdW t + Vtdt. 



< 




(6) 



dX t = TlX t dt + ZdW t + rjtdt, 
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(1) Principal Components analysis (PCA) on the returns (they are computed 

, F-(t+l day,T)-F e (f,TK 
uy Fe(t,T) r 

(2) Estimation of the parameters rf (i = 1, 2) with a nonlinear regression 

T-t T—t 

between 1, e , ^-^e and the motions X deduced from the PCA. 

(3) Reconstruction of the differences dXf at each time step with linear re- 
gression between the volatility functions vector and the returns vector. 

(4) Linear regressions between (AX t ) and (X t ) to determine the matrices IT 
and ££*. In the matrix II, we keep significant elements with a regression 
subset selection according to the Bayesian Information Criterium (BIC) , 
see Miller Q. 

(5) Computation of the function (r} t ) t . 

Only the last step needs to be detailed. We mention first that we have set 
Xq = in the previous steps, which is equivalent to shift the (unknown) value 
of (r] t )f 

(i}t)t is an additional drift parameter. Usually, it can be identified only over a 
very long time data set, which is impossible in practice. 

To overcome this undertermination problem, we propose to adjust it so that 
expected forward prices will fit the initial curve (for each energy): 



This is often expected by practitioners. Note that this adjustment implies 
a vanishing market risk premium 7r e (0,T) but any alternative choice would 
be possible. The model is now rewritten in a more simple manner. We set 
Xt = Xt + 9 t where (9 t )t solves 9 t = f (H9 s + ^s)ds, which is given by 
9 t = e m Jq e~ sU r] s ds. Thus, it is equivalent to determine r] or 9. Equations 
(6) consequently become 



E (F e (t, T)) =F e (0,T), 



V(t,T). 



< 



dF3(t,T) 
F9(t,T) 

&F c (t,T) 
F c (t,T) 



= a9(T-t)(dX t + 9'tdt), 
= a c (T -t)(dX t + 9' t dt), 

Tix t dt + mw t} X = 



(7) 



dX t = 



0. 



In the Appendix, we prove that 




In addition, a closed formula for Eg =0 (fwof^) is given (see (10)) which only 
depends on the already estimated parameters. In order to have the left hand 
side equal to 1, one should have 





dt 
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for any (t, T) and for each energy, with a suitable choice of 9. In our tests on 
market data, the perfect fit is achieved with less that 0.1% error (see Figures 
7). In that case, 9 is obtained by a least squares optimisation routine applied 
to the difference of the two sides of (9). 

Note that when IT is zero, the forward prices are already centred on the initial 
forward curve if 9 = 0: no adjustment of 9 is needed. X are just non standard 
Brownian motions and the model comes down to the one given by Brooks. 



3 Numerical simulations 

The historical data comes from the ICE (InterContinentalExchange) Market. 
The quotations begin in September 8th, 2003 and end in april 5th, 2007. 9 
contracts for natural gas prices and 15 for crude oil are quoted for different 
maturity months. 

The Figure 2 represents the prices of the contract delivering one unit of gas 
or one unit of oil the next month (on the right hand side) and nine months 
later (on the left hand side). Even if the gas time series is seasonal, we notice 
that natural gas and crude oil prices have the same long term tendency. 



One month contract for gas and oil Nine month contract for gas and oil 




2004 2005 2006 2007 2004 2005 2006 2007 

quotation dates quotation dates 



Figure 2. Forward contracts delivering one unit of natural gas or crude oil in one 
month or nine months (p/th for natural gas and $/bbl for crude oil). 
Despite the seasonality of gas, natural gas and crude oil prices seem to have a 
common long term tendency. 

To avoid the seasonality effect, we can plot the gas and oil prices for a fixed 
delivering date (fixed maturity), see Figure 3. This graph illustrates the long 
term tendency of natural gas and crude oil. 

The first three steps of calibration allow us to obtain the following parameters 
(for N 9 = N c = 3): 
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Contract delivering in jan 05 



Contract delivering in jan 05 




may Jul. sept. nov. 1 1 1 r~ 

40 50 60 70 

quotation dates (Irom avr 04 to dec 04) 



Contract delivering in sept 05 Contract delivering in sept 05 




jan. mar. may Jul. I 1 1 1 r~ 

26 2S 30 32 34 

quotation dates (from dec 04 to aug 05) 



Figure 3. Contracts delivering in January and September 2005 (p/th for natural gas 
and $/bbl for crude oil). 

Natural gas and crude oil prices have a common long term tendency. 



r! 


r 9 




'2 


0.736 year 


0.086 year 


3.761 years 


0.138 year 



The parameters [j^) e convey the characteristic time of a short term shock. 
The impact of a shock on crude oil prices will last during a longer period than 
a shock on natural gas prices (i.e. 3.761 years > 0.736 year). 

Using Brooks' model with a simple correlation between Brownian motions to 
link energies, we can obtain such prices' simulation as in Figure 4. 

From the market data and the previous values of rf , we obtain the motions 
X t (see Figure 5). 

Figure 6 shows an example of stationary process resulting from one of the long 
term relations binding the motions {X t ) t . 

Cointegration tests (Phillips-Ouliaris) allow us to reject the null hypothesis of 
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Simulation of the contract delivering in dec 07 




Simulation of the contract delivering in dec 07 

55 60 65 70 




quotation dates (from mar 07 to nov 07) 



50 55 



Figure 4. Simulation for the contract delivering in December 2007 using Brooks' 
model. 

Natural gas and crude oil prices have no common long term tendency. 



Motions Xt for gas Motions Xt for oil 




2004 2005 2006 2007 2004 2005 2006 2007 

quotation dates quotation dates 



Figure 5. Motions (Xt)t for natural gas and crude oil. 

According to the Phillipss-Ouliaris cointegration tests, the motions (Xfjt are coin- 
tegrated. 



non-cointegration. This justifies the use a Vectorial Error Correction Model 
to model the motions X t . With significant linear regressions (according to the 
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Stationary process of a long term relation 




2004 2005 2006 2007 

quotation dales 



Figure 6. Stationary process from one of the long term relations binding the motions 
X t . 

The long term relation presented here is X^ — 0.97 Xf = ut where (ut) is a stationnay 
process. 

BIC), we find the following matrices for IT and EE*: 



IT = 



EE* = 





-0.017 


0.019 


) 











0.005 0.009 -0.027 


-0.162 











-0.012 


0.174 











-0.009 


-0.030 











0.015 


0.008 0.046 






V 





-0.017 


0.019 -0.052 ) 






/ 


0.00158 


-0.00323 0.00386 


-0.00001 0.00006 


0.00003 N 




-0.00323 


0.00812 -0.00958 - 


-0.00007 0.00007 


-0.00007 




0.00386 


-0.00958 0.01740 


0.00006 -0.00003 


0.00011 




-0.00001 


-0.00007 0.00006 


0.00045 -0.00052 


0.00010 




0.00006 


0.00007 -0.00003 - 


-0.00052 0.00096 


-0.00011 


V 


0.00003 


-0.00007 0.00011 


0.00010 -0.00011 


0.00015 j 



All these parameters allow us to compute E^ =0 (j^o?)) * n orc ^ er to centre the 
simulations on the initial forward curve. Figure 7 shows this result for fixed 
(T — t). We now show some examples of price simulations using this model. To 
simulate one scenario, we use an Euler scheme. The prices are diffused from 
the initial forward curves of March 15th, 2007. We draw the simulations for a 
fixed maturity so as to favour the long term relation expressed by the model, 
over the seasonality of gas. Three simulations are shown on Figure 8. These 
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E (F(t,T)/F(0,T)) for theta(t)=0 
Gas 



E (F(t,T)/F(0,T)) for theta(t)=0 
Oil 



■ Month 1 
Month 2 
Month 3 
Month 4 
Month 5 
Month 6 
Month 7 
Month 3 
Month 9 



■ Month 1 
Month 2 
Month 3 
Month 4 
Month 5 

■ Month 6 
Month 7 
Month 8 
Month 9 



Jul Sep 
quotation dales 



E (F(t,T)/F(0,T)) for theta(t) adjusted 
Gas 



■ Month 1 
Month 2 
Month 3 
Month 4 
Month 5 
Month 6 
Month 7 
Month 8 
Month 9 



E (F(t,T)/F(0,T)) for theta(t) adjusted 



■ Month 1 
Month 2 
Month 3 
Month 4 
Month 5 

■ Month 6 
Month 7 
Month 8 
Month 9 



Figure 7. Thanks to the parameter 8t, prices' simulations are centred on the initial 
forward curve. 

The ratio of expectations must be 1. In the first case (it, 9% = 0), we notice a 
maximum difference with 1 of 5.6% for gas and 2.1% for oil. In the second case (9t 
is now adjusted), the maximum deviation from 1 is 0.08% for gas and 0.03% for 
oil. 

simulations fit well the possible scenarios given by market data (see Figure 3). 



Conclusion 



In this work, we propose a continuous time model for natural gas and crude oil 
future markets, which brings together two key features of these energy markets. 
On the one hand, it conveys a long term dependence between gas and oil, as 
it has already been supported by several previous econometric studies. This 
dependence is achieved through a cointegration of motions driving the prices. 
On the other hand, the model is free of arbitrage and coincides with usual 
forward models under the risk neutral probabilities. As a consequence, this 
model can be used for Value at Risk computations (or other risk management 
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Simulation 1 of the contract delivering in dec 07 



Simulation 1 of the contract delivering in dec 07 




40 50 60 70 




quotation dates (from mar 07 to nov 07) 



20 25 30 35 40 45 50 



Simulation 2 of the contract delivering in dec 07 




Simulation 2 of the contract delivering in dec 07 

55 60 65 70 




quotation dates (from mar 07 to nov 07) 



Simulation 3 of the contract delivering in dec 07 




Simulation 3 of the contract delivering in dec 07 

GO 70 SO 90 100 




0.87 











4b 





crude 



50 60 70 SO 3D 



quotation dates (from mar 07 to nov 07) 



Figure 8. Three examples of simulations for contracts delivering in one month. 
The model conveys the long term relation between natural gas and crude oil con- 
tracts. 

measurements) as well as for option pricing. 

Based on market data related to future contracts, the calibration gives very 
good results. It would be presumably improved by incorporating derivatives 
prices in the calibration set. 

We have presented a version of the model without seasonality effect (especially 
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regarding gas). We mention finally that it is easy to incorporate seasonality 
in the gas volatility. 
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Appendix 



We first recall the exponential form of the forward price: 
F e {t, T) = F e {0, T) exp a e {T - s)9' s ds 



exp 



o 



* a*(T - s)dX s -\(<j e {T - .)X) t ) . 



The above quadratic variation equals / * \a e (T — s)E| 2 ds. Since X does not 
depend on 9, the assertion (8) is proven. 



Let us prove the following formula 



E 



6=0 



'F e (t,ry 



■ exp 



Jo {l ll^^ ~ M ) ne( " _S)IIdMS 



(10) 



+ a e (T - s)££* J* e (M " s)n *n*(a e (T - u))*duj ds 



Set I = j* * a e (T - s)dX s ; then, one has 



E- 



P (F*(t,T) 
-° \F e (0 



^) = E P (exp(J)) exp (-1 jf K(T - S )S| 2 d S ) . (11) 



Our methodology to compute E p (exp (J)) consists in showing that I is a cen- 
tred Gaussian variable with a variance v to calculate; it follows that E p (exp (J)) 
exp(t> 2 /2). The process X given in Equation (7) is a generalized Ornstein- 
Uhlenbeck process. Its solution is explicit and is given by: 



X t = C e (t - s)n EdW s . 
Jo 

Then, using again Equation (7), we obtain 

/= /* a e (T -s)UX 9 ds + f a e (T - s)mW s 
Jo Jo 

= I a e {T-s)U f e {s " u)u mW u ds + [ a e (T - s)mW s 
Jo Jo Jo 

= J* J* o e (T - u)Ue iu - s)u mu + o e (T - s)E 



(12) 



dW« 



Thus, J is a Wiener stochastic integral, whose is centred and whose variance 
equals / * |/* a e (T - ^Ue^-^du + a e (T - s)s| 2 ds. Changing in (11) and 
making simplifications lead to (10). 
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One month contract for gas and crude 
Simulation 1 




One month contract for gas and crude 
Simulation 2 




One month contract for gas and crude 
Simulation 3 




One month contract for gas and crude 
Simulation 4 
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